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A re-evaluation of time-averaged accretion rates at DBZ-type white dwarfs points 
to historical, time-averaged rates significantly higher than the currently observed 
episodes at their DAZ counterparts. The difference between the ongoing, instanta- 
neous accretion rates witnessed at DAZ white dwarfs, which often exceed lO^gs"^, 
and those inferred over the past 10^ — 10^ yr for the DBZ stars can be a few orders of 
magnitude, and therefore must result from high-rate episodes of tens to hundreds of 
years so they remain undetected to date. This paper explores the likelihood that such 
brief, intense accretion episodes of gas-phase material can account for existing data. 
For reasonable assumptions about the circumstellar gas, accretion rates approaching 
or exceeding lO^^gs"^ are possible, similar to rates observed in quiescent cataclysmic 
variables, and potentially detectable with future x-ray missions or wide-field monitor- 
ing facilities. Gaseous debris that is prone to such rapid accretion may be abundant 
immediately following a tidal disruption event via collisions and sublimation, or if 
additional bodies impinge upon an extant disk. Particulate disk matter accretes at or 
near the Poynting- Robertson drag rate for long periods between gas-producing events, 
consistent with rates inferred for dusty DAZ white dwarfs. In this picture, warm DAZ 
stars without infrared excesses have rates consistent with accretion from particulate 
disks that remain undetected. This overall picture has implications for quasi-steady 
state models of accretion and the derived chemical composition of asteroidal debris in 
DBZ white dwarfs. 
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1 INTRODUCTION 

Atmospheric metals should not be present in isolated 
white dwarf stars with effective temperatures below roughly 
25 000 K. At this stage of their evolution, radiative forces 



become insignificant (jChaver et al.|[l995l ). and gravitational 
settling pulls elements heavier than helium into the stellar 
interior on timescales a few to several orders of mag nitude 
shorter than their cooling ages (|Fontaine et al.|[l979l ). Nev- 
ertheless, white dwarfs in this temperature range exhibiting 
atmospheric metal absorption have been known for nearly a 
century, and the pro totype is still the nearest single white 
dwarf known, vMa 2 (|van Maanen|[l917l ). At 4.4 pc, this star 
exemplifies a few characteristics of the metal pollution phe- 
nomenon: it is a single star and thus not accreting from 
a stellar or substellar companion ijFarihi et "ai] |2008al ). its 
location w ithin the Local Bub ble and helium-dominated at- 
mosphere (|Dufour et al.llioOTi ') precludes accretion from in- 
terstellar material, and its age of 3 Gyr is at odds with the 



persi stence of metals that sink within a few Myr l|Koesterl 
[2003). Combined with the above facts, the abundance pat- 
terns in stars like vMa 2 have long suggested rece nt accretion 
from material such as that found in asteroids (|Sion et al.l 
l1990; Graham ct al. 1990). 

Indeed, evidence is now strong that metal-polluted 
white dwarfs accrete from planetary debris. Spitzer stud- 
ies have revealed circumstellar disks at a significant frac- 
tion of these stars, where the temperature profile and in- 
ferred geometry of the dust is consistent with material com- 
pletely contained within the Roche limit for km-sized or 
larger solid bodies, with innermost disk temperatures that 
shoul d rapidly sublimate solids, and drive material onto the 
star jXu fc Jural I2OIII: iFarihi et al.ll2009l: ljura et al.l 12007 



Hipp el et all l20o'^ : iGansicke et all I2OO6I : iReach et all 
2005). The favored model for the origin of thi s material is 
the tidal disruption of a large asteroid analog (|jura|[20o3 '). 
perturbed into a high eccentricity by an unseen body such 
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as a maior 



planet 



iDebes fc SigurdssorJ 



Debes et al] I2OI2I : iBonsor et all I2OIII : 



20021 ). All available mid-infrared spec- 



trosc opy reveals disk matter that is silica te-rich and carbon- 
poor (|jura et al.ll2009l : lReach et al.ll2009l ). while stellar spec- 
troscopy reveals atmosphe ric pollution by refractory-ric h 
and volatile-poor material (|Farihi et al.ll20l"0al : |jurall2006l ). 
Critically, high-resolution optical spectroscopy of the disk- 
contami nated stars r eveals distinctly terrestr ial-like abun- 
dances l|Klein et alj I2OI0I : IZuckerman et al~ I [2007). Thus 
metal-enriched white dwarfs are astrophysical traps for exo- 
planetary debris, acting as detectors that can yield the bulk 
composition of planetary building blocks that orbit interme- 
diate mass stars. 

This paper re-examines in detail the inferred accre- 
tion rates at both hydrogen and helium atmosphere, cool 
white dwarfs with trace metals, referred to here as DAZ and 
DBZ stars resp ectively, based on st andard spectral types. 
As described in iGirven et al.l l|2012l ') , an improved calcula- 
tion method is used to place both classes on equal footing, 
and in doing so, the derived accretion rates, averaged over 
the recent history for the DBZ stars, are found to be signif- 
icantly higher than the current rates observed at the DAZ 
stars. This work explores these issues in detail and develops 
a model to address the observations. In §2 the accretion rate 
calculations are described and examined, and §3 presents a 
scenario that can account for these data. Implications for 
future observations and the nature of the destroyed and ac- 
creted planetary bodies are discussed in §3 and §4. 



2 INFERRED ACCRETION RATES FOR DAZ 
AND DBZ STARS 

2.1 Steady State Metal Accretion 

If a star is in a steady state between accretion and diffu- 
sion, the total mass accretion rate for heavy elements can 
be expressed as (|Dupuis et al.lll993l ) 



(1) 



where Xi is the mass fraction, the diffusion (i.e., sinking) 
timescale for the i^^ heavy element, and Mcvz is the mass 
of the stellar convection zone or outer atmosphere (defined 
at the base of the convection zone or at Rosseland optical 
depth 5, whichever is deeper; [Koester 2009). The numerator 
on the right-hand represents the mass of a the i**" element 
in the convective or outer layers of the star. 

All known metal-polluted white dwarfs exhibit the Ca 11 
K line in optical spectroscopy owing to the strength of this 
atomic transition at this range of effective temperatures (for 
this reason it is the strongest absorption feature in the Sun). 
Accretion rate calculations have been naturally tied to this 
element, and one can derive an accurate infall rate for cal- 
cium, and extrapolate to the total mass accretion rate with 
an appropriate correction factor, as follows 



A 



T"Ca 



(2) 



Prior calculations assumed the infalling material 
was exactly solar (i.e., including hydrogen and helium; 



iKoester fc WilkenlfioO^. or 1% sol ar by mass (i.e., no hy- 
drogen or helium; Jura et al]l2007l ): the latter yields A = 
1/109.1 or 1/42.7 depending on whether the extrapolation 
is made from a calci um abundance measured as [Ca/H] or 
[Ca/He] respectively (lLoddersll2003l ). 



2.2 Instantaneous vs. Historical Rates 

For Teff > 11 000 K DAZ white dwarfs. Equation [T] can be 
applied with confidence because these stars possess tiny con- 
vection zones or none at all, and have commensurately short 
metal sinking timescales (on the order of days; JCocstcr 200^ 
that imply a steady state is virtually certain. Rates com- 
puted in this way for warm DAZ stars are ongoing and in- 
stantaneous measures of the infall of a given heavy element. 
In contrast, for the coolest DAZ stars and the DBZ stars in 
general, the assumption of a steady state is uncertain owing 
to larger diffusio n timescales from cor respondingly deeper 
convection zones ijPaauette et al.l[l986l ). The DBZ stars can 
retain atmospheric metals for timescales up to 10^ — 10^ yr, 
and while they cannot yield instantaneous accretion rates, 
their sizable convection zones contain the integrated signa- 
tures of accretion; a substantial mass of metals acquired 
within the last few sinking timescales. Because of this, DBZ 
stars provide a historical record of accretion and, impor- 
tantly, a minimum mass for the destroyed parent body or 
bodies (jKlein et al.ll2010l : iFarihi et al.ll2010al : Ijurall200d l. 

In ord er to best c o mpar e these two classes of pol- 
luted stars, iFarihi et al.l (|2009l ) calculated historical accre- 
tion rates for DBZ stars; that is, time-averaged over a sin- 
gle diffusion timescale, assuming a steady state. The like- 
lihood of an infrared excess due to circumstellar dust at 
DAZ stars is correlated with ongoing, high-rate accretion, 
and thus a similar correlation in the DBZ stars was natu- 
rally sought. While the fraction of stars with infrared excess 
differs between the DAZ and DBZ white dwar fs with accre- 
tion r ates above a benchmark of 3 x 10^ gs ~^ (iGirven et al.l 
l2012l : IXu fc Jural I2OI1I : IFarihi et al.ll2010bl . |2009| ). the rele^ 
vant point for this study is that until recently, in terms of 
metal pollution alone and based on the assumption of so- 
lar calcium relative to other heavy elements, the two stellar 
classes appeared to have broadly similar, inferred metal ac- 
cretion rates. 

In what follows, the terms instantaneous and ongoing 
are used to refer to accretion rates for warm DAZ stars where 
a steady state is highly probable, while historical and time- 
averaged refer to the rates derived for DBZ stars and the 
coolest DAZ stars, where an integrated approach is neces- 
sary. 



2.3 New Accretion Rate Calculations 

Assuming 0.01 solar composition by mass, extrapolated from 
calcium abundance, for both DAZ and DBZ white dwarfs is 
not ideal, but understandable based on exist ing data just 
a coup le years prior. However, as discussed in IGirven et al.l 
(j2012l) . a more accurate estimation of accretion rates is pos- 
sible based on recent work. Briefly, there are now eight 
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Figure 1. Total metal accretion rates onto 78 contaminated stars 
observed with Spitzer IRAC. The open and filled symbols respec- 
tively denote DAZ and DBZ-type white dwarfs, while the stars 
and circles respectively distinguish those white dwarfs with and 
without infrared excesses due to circumstellar dust. Among these 
targets alone there are 8 DBZ stars with accretion rates above 
the maximum witnessed to be ongoing for any known DAZ white 
dwarf. Cooling ages are calculat ed baaed on white dw arf atmo- 
spheric and evolutionary models llFontaine et al.ll2001^ . 



stars with measured O, Mg, Si, Ca, and abundances, 
with some of these s tars exhibiting further elements such as 
Al, Ti, C r, and Ni ('Gansicke et al."2012'; 'Klein et al."201l|; 
Melis et al. 2 011: Farihi ct al. 2011a; Zuckcrman ct al. 2gip; 



Vennes et al. I bold iDufour et all I2OI0I : iKlein et all I201C 



These data demonstrate that the accreted matter in polluted 
white dwarfs is broadly simila r to rocky material of th e in- 
ner Solar System. Specificallv. IZuckerman et al.l (|2010l l find 
that calcium represents, on average, close to 1:60 by mass, of 
all the accreted heavy elements for well-studie d stars, while 
this same ratio is 1:62.5 for the bulk Earth (|Allegre et al.l 
1 19951 ). 

These results suggest an improved metal accretion rate 
for dozens of stars where only calcium is detected can 
be made assuming it is 1.6% {A = 1/62.5) of the to- 
tal mass. Figure [T] plots accretion rates calculated using 
Equation [2] fo r 40 DAZ arid 38 DBZ stars observed with 
Spitze r IRAQ jGirven et al.|[20]l : iFarihi et al.ll2010bl . |2009| . 
l2008bl : lJuraet al.l |2007| ). Calcium abundances and stel- 
lar parameters come fr om the most recent studies (e.g., 
iKoester fc Wilkenll2006l ). with si nking timescales and con- 
vective envelope masses given in iKo ester (2009). Compar- 
ing this plot with those published previously, it is clear there 
are an increasing number of outlying DBZ-type stars with 
inferred rates one to a few orders of magnitude above that 
observed at any DAZ star. Figure [2] re-plots these accre- 
tion rates as histograms for both atmospheric types (ig- 
noring dust emission), clearly delineating a subgroup of 
high values that occur exclusively for DBZ stars. While 
GALEX 1931-f0117 is currently the most extreme exam- 
ple of a DAZ sta r, with an instant aneous accretion rate of 
2 X 10^ gs"^ (|Gansicke et al.ll2012l '). its extreme DBZ-type 



^ These elements co mbined represent over 95% of the bulk Earth 
llAllegre et al.lll995h 



Figure 2. Histogram of the rates plotted in Figure [T] for each 
atmospheric class of polluted white dwarf. The population of rates 
log{dMz /dt (gs)~^) > 9.5 occur only among the DBZ-type stars. 



counterpart is HE 0446— 2531 which has a time-averaged 
rate two orders of magnitude higher at 3 x 10^^ g . Figures 
[T] and [2] plot a large subset of all nearby (pre-SDSS) metal- 
enriched white dwarfs, but there are add itional examples o f 
high rate DBZ stars (e.g., Figure 11 of IXu fc Jural 120111 ), 
and at least 10 cases where the historical accretion rates are 
higher than any inferred to be ongoing at DAZ stars. 

Table [1] compares accretion rate calculations for the 
eight stars with robust detections and abundances of O, 
Mg, Si, Ca, and Fe: in the second column using Equation 
[2] and assuming calcium is 1.6% of the total accreted metal 
mass (as described above), and in the third column from 
the actual detected elements using Equation [T] While a few 
of these stars exhibit additional elements beyond the five in 
common, the third column values represent a minimum, as 
the detection of further elements (e.g., Al, Ni) may serve to 
increase the calculated rates. However, because most of the 
major terrestrial elements are detected in these stars, it is 
possible but unlikely their total accretion rates will change 
substantially. Comparing the inferred rates from both meth- 
ods, there is agreement to within ±0.2 dex for five of the 
eight stars and overall, while three are over-predicted. These 
over-predictions may diminish as further elements are de- 
tected in these stars, but it is more likely that the (modest) 
disagreement in the calculated accretion rates are due to 
real metal-to-metal abundance variations among the parent 
bodies now falling onto the stars as debris (|Zuckerman et aP 
I2OIII : iKlein et al.ll201Gl ) . 



2.4 Potential Biases 

It is worth noting that GALEX 1931 is the only known DAZ 
star with optical detections of all the necessary metals for 
the Table[T]comparison. This is because, all else being equal, 
the atmospheric opaci ty of hydrogen is signifi cantly higher 
than that of helium (jZuckerman et akl |2003| ) . Because of 
this, there is a bias against the optical detection of DAZ 
star s with modest to low am ounts of pollution (see Figure 
1 of IKoester fc WilkenllioO^ for an excellent illustration of 
this fact). Within the DA cooling sequence, the combination 
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of atmospheric opacity and Ca ll ionization fraction permits 
lower metal abundances and accretion rates to be detected 
at cooler stars. But between the two classes, the detection 
of metal lines in DA stars is more difficult, favoring those 
metal-rich stars exhibiting relatively stron g lines, and thu s 
high abundances and accretion rates (|Koester et al.|[2005h . 
Because the bias works in the opposite direction, the lack of 
high rates in the DAZ population is real. 

Another possibility is that more vigorous accretion 
takes place at DBZ stars. However, it is difficult to con- 
jure any scenario that physically favors higher rate accre- 
tion of planetary debris at stars with helium atmospheres. 
There is neither empirical nor theoretical reasoning that sug- 
gests the existence of DB stars is related to anything be- 
yond (single) stellar evolution: a very late thermal pulse, 
immediately following the asymptotic giant phase, can en- 
rich t he outer layers of a star with hydrogen-deficient ma- 
terial l|Herwig et al.|[l999|). event ually producing a helium- 
dominated remnant (|Werner fc Ranch 1994). Other expla- 
nations along these lines would include that DB stars tend 
to have planetary system architectures more efficient at as- 
teroid perturbation, or parent bodies highly enriched in cal- 
cium, but such solutions are contrived and therefore dis- 
counted. 

Lastly, white dwarf atmospheric models do not have 
substantial uncertainty in the sizes of the convective or outer 
envelopes, nor in the gravitational s ettling behavior at the 
bottom of these zones ([Koes tcr 2009). However, it is feasible 
that there are gaps in the theoretical understanding of stel- 
lar atmospheres that may contribute to the inferred trends, 
e.g., if the settling times in Equation [1] are systematically 
under-predicted for DB stars. This possibility is outside of 
the scope of this paper and not discussed further. 

The following sections proceed on the assumption that a 
planetary system forms without any foreknowledge of the ul- 
timate stellar remnant, and that the planetary bodies within 
depend only on the conditions at formation and the subse- 
quent influence of stellar evolution and mass loss on those 
bodies (i.e., only stellar radiation and dynamical changes). 



3 A HISTORY OF HIGH ACCRETION RATES 

The inescapable conclusion of the data shown in Figure [5] is 
that there must be relatively short-lived episodes of accre- 
tion at metal-polluted white dwarfs; brief enough not to have 
been observed as ongoing, and at rates substantially higher 
than those witnessed to date at DAZ stars. This is obvi- 
ously true for the six stars with time-averaged rates above 
lO^^gs"^, but it also appears likely for the greater DBZ 
population. Comparing the rates between the DAZ and DBZ 
groups plotted in the figures, the highest accretion rates dif- 
fer by a factor of 160, the average of the top six rates diff'er 
by a factor of 40, and the average of all the rates differ by 
a factor of 25. Thus the phenomenon is represented by the 
entire class of helium atmosphere stars, and accentuated by 
a handful of outstanding examples. 

The rates inferred for the DBZ stars are the mass of 
accreted metals over a single diffusion timescale, divided by 
that time scale; around 10^ yr at 17 000 K and 10® yr near 
12 000 K (iKoesteijlioOol ). However, in cases where circum- 
stellar dust is not detected at a given DBZ white dwarf, it 



Table 1. Comparison of Accretion Rate Estimates For Stars with 
Detected O, Mg, Si, Ca, Fe 



Star 






G241-6 


9.8 


9.3 


GALEX 1931 


9.2 


9.2 


GD 40 


10.1 


9.4 


GD61 


9.0 


8.8 


HS 2253-1-8023 


10.2 


10.1 


PG 1015+161 


8.5 


8.2 


SDSS0738 


10.3 


10.3 


SDSS 1228 


9.0 


8.8 


Log Average 


8.66 


8.43 



" Assuming Ca is 0.016 of the total accreted mass. 

^ Fr om measured a bunda nc es of at le a st O, Mg, Si, Ca, 

Ffe I Gansicke et al.l |2012| : iKIein et al.1 l201ll: iFarihi et al 

'2011a : lMelisetal.1 |2010|: IVennes et al.l |2010| : iKlein et a" 
l201(rlDufour et al.ll2010l : IZuckerman et al.ll2010l) . 



Note. The overall differences in the inferred accretion rates be- 
tween the two columns are minor, and will not affect the findings 
for DBZ versus DAZ properties; see §2.3 for calculation details. 

is possible that accretion halted within the past one to sev- 
eral diffusion timescales. In this sense, the rates inferred for 
the DBZ stars are minimum, time-averaged rates, as atmo- 
spheric metal abundances will exponentially diminish, but 
remain potentially detectable for a few to several sinking 
timescales after accretion has ended. For simplicity, the fol- 
lowing assumes all observed metals are accreted over a sin- 
gle diffusion timescale; if accretion is no longer ongoing, then 
prior accretion must have occurred at even higher rates than 
those inferred below. 



3.1 High and Low State Model 

For a toy model with two distinct and constant accretion 
rates - a high state and a low state - taking place over a 
single diffusion timescale for a typical DBZ star, the ob- 
served mass of metals in the stellar convection zone Mz can 
be written as: 



MzT = Mhth + Miti 



(3) 



where t^/r is the fraction of time the star is accreting at 
the high rate Mh, similarly ti/r is the fraction of time spent 
accreting at the low rate Mi, and where r is the diffusion 
timescale. The right hand side of Equation [3] simply breaks 
the total accreted mass into that delivered in two distinct 
states, but one can easily extend this to a continuum of 
rates and corresponding intervals; this does not change the 
conclusions reached here. As shown in §3.2, since th <^ti, it 
follows that ti ^ T and thus t^/ti is the fraction of the total 
time spent accreting in the high state. This allows for the 
possibility that Mh occurs in a succession of episodes, not 
just one per DBZ diffusion timescale. The mass represented 
by both sides of Equation [3] is typically 10^^ g for a given 
DBZ star, but this can v ary by up to two orders of magnitude 
(|Girven et a l.| [2012 : Farihi et al'lliblOa ). 

Identifying r and Mz as the diffusion timescale and 
time-averaged accretion rate for DBZ-type stars, and equat- 
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Table 2. Accretion Rates for Figure 3 





DBZ stars 


DAZ stars 




Line 


log(M,[gs-l]) 


log(A-/ags-i]) 


Description 


Dash 


11.5 


9.3 


Max M 


Long Dash 


10.8 


9.2 


Top 6 M 


Solid 


10.1 


8.7 


Avg M 



Note. There are six DBZ stars with inferred accretion rates above 
10^'gs^l (see FigureEJ. 

ing Ml to that measured for the DAZ white dwarfs, Equation 
[3] permits the assessment of {th,Mh) values sufficient to ac- 
count for the inferred differences among the instantaneous 
and historical accretion rates. Figure [3] plots three represen- 
tative solutions for r = lO'' yr, with parameters described 
in Table (2] From top to bottom are plotted the resulting 
(thjMh) for pairs of Mz and M; representing 1) the single 
highest accretion rates in each class, 2) the average of the 
top six rates in each class, and 3) the average among all stars 
in each class from Figure [T] 



3.2 Allowed Timescales for High Rate Episodes 

Because ongoing accretion rates as high as those inferred 
for the DBZ have not been witnessed to date, strict limits 
can be placed on the timescales for such episodes, th- For 
r = 10® yr and a high-rate accretion timescale of th = 10^ yr 
(O.llti), the probability that at least one of 46 known DAZ 
stars will be accreting at the high rate is above 99%, while 
for th = 10* yr (O.Olt;) this same probability drops to 37%. 
Based on this argument alone, one can conclude th "^ti, and 
the timescale for any high-rate episodes must be no more 
than 10* yr. Moreover, a few of the most highly polluted 
DBZ stars have r ~ 10^ yr (e.g., GD 362, Ton 345), which 
supports even shorter timescales for high-rate spikes. 

Notably, DA stars without detected metals constrain 
th m ore strongly. The high- resol ution spectroscopic s tud- 
ies of lZuckerman et al.l (|2003l ') and lKoester et al.l ^200^ col- 
lectively surveyed approximately 530 DA white dwarfs for 
photospheric metals, where stars more highly polluted than 
GALEX 1931 would have been readily detected. It is also 
likely that extreme metal abundances in DA white dwarfs 
are identifiable in the low-resolution spectra of the Sloan 
Digital Sky Survey; SDSS 1228 is one of the most metal- 
rich DAZ stars and its discovery spectrum e xhibits a clear 
Mg 11 absorption line at <; = 16.2 AB mag (|Gansicke et al.l 
I2OO6I ). While this star is relatively bright among the few 
thousand spectroscopically confirmed DA white dwarfs iden- 
tified in SDSS DR7, there are nearly 470 DA stars with 
g < 17.0 AB mag where very high met al abundances (an d 
thus accretion rates) were not detected (|Girven et al.|[201ll ). 
Based on these surveys, it is likely that th < yr. 



3.3 Rapid Accretion of Gaseous Debris 



Using the a prescription I King et al.l |2007| ) for completely 
gaseous accretion disks, 'jura' (200^) estimates that a disk 
composed purely of metallic gas will dissipate within 10^ yr, 
for a — 0.001. While this timescale is longer than the allowed 
range of th in Figure [31 below it is shown that a reasonable 



10 




Detection Probability: 10% I 50% I 100% I 



Highest Observed Ongoing dM/dt 



log[t^ (yr)] 

Figure 3. Necessary (th,Mh) that satisfy Equation [3] for r = 
10^ yr, and Table |2] values for and M;. It should be noted 
that there is little difference in the plotted lines and those that 
result from Mi = 0. Also shown are resulting benchmark prob- 
abilities of detecting high states of marked durations, given the 
sample of 534 DA stars surveyed at high resolution with Keck 
and VLT llKoester et al.ll2005l : IZuckerman et al.ll2003l) . Episodes 
of high rate accretion lasting tens to hundreds of years are al- 
lowed, while if the timescales approaches or exceeds 10^ yr then 
the likelihood of detection at DAZ stars increases towards unity. 
This result is corroborat ed by the lack of ex tremely polluted DAZ 
stars in the SDSS DR7 dCirven et al.ll201ll) . 



choice of a can produce both th and Mh values necessary to 
account for the accretion rate data. 

The accretion rate of a completely gaseous disk can be 
expressed as 



M = — 

t^ 



(4) 



where m is the mass of the disk and t^ is the viscous 
timescale at a distance d from the central star 



(5) 



The viscosity v can be rewritten following the a prescription 
IPringleHigSlI ) and Equation [5] becomes 



(6) 



where p is the orbital period and h is the scale height of the 
disk. For an ideal, isothermal gas o rbiting a star of mass M 
at a distance d, the scale height is (|Pringlelll98"ll '): 



h = 



kTd? 
GM/i 



(7) 



where T is the temperature of gas with mean molecular 
weight fi. 

Metallic, circumstellar gas has been modeled at temper- 
atures around 6000 K in the few, well-studied stars with cal- 
cium triplet emis s ion li n es, using three ind e pendent means 
(Hartman n et al.l 1201 ll : iMehs et al] l2010l : iGansicke et all 
2006). Gas composed of singly ionized O, Mg, Si, Fe will 
have a mean particle mass fj, = 2.5 x 10"^'^ g, and if this or- 
bits at 0.5 Rq about a 0.6 Mq white dwarf, then d/h ~ 260. 
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Then the viscous timescale at this distance becomes ~ 
(10/Q)yr and in the range 25 to lOOyr for a = 0.1 — 0.4; 
typically required to account for observational data in a wide 
variety of fully ionized disks (|King et al.ll2007l ). This viscous 
timescale is consistent with observations in two ways: 1) it 
is sufficiently short as to remain currently undetected at 
DAZ stars, and 2) for the largest, acc reted metal masses 
in DBZ star conve ction zones flO^** g: Icirven et al.l [201^ : 
iFarihi et al.l l2010al ). it implies rates up to lO^^gs"^ (and 
higher for smaller d). 

Metal-rich, pure gas disks at such temperatures will be 
highly or fully (singly) ionized as their overall ionization 
potential will be lower than similar accretion disks of so- 
lar abundance, hydrogen-dominated materiaffj. Thus, a disk 
of gaseous metals wil l be conductive, s t able, and the stan- 
dard a should hold (|Stone et all I2OO0I : [ Gammie fc Menoul 
Il998l ). Accretion of gaseous debris via viscous dissipation 
with a — 0.1 — 0.4 thus nicely accounts for the data in 
Figure [2I yielding the correct range of rates and timescales 
dictated by Figure [S] and is hence plausible. 

Two recent studies have also shown that short bursts 
of high rate accr etion are possib le from disrupted asteroids 
at white dwarfs. iRafikovl l|2011bl ) demonstrates that, under 
certain conditions, gas may efficiently couple with solids 
and lead to the inward drift of particulate disk matter at 
rates greatly exceeding that of Poynting-Robertson forces. 
In this scenario, a spike o f maximal accretion occurs at the 
very end of a disk episode (Mctz ger et al.|[2012l ). In contrast, 
iBear fc Soken (2012) show that a Nova-like, accretion-driven 
outburst may occur during the initial phases of an asteroid 
disruption and disk formation. The standard a disk model 
used here applies to any epoch where sufficient gas is pro- 
duced, and overlaps with these other models under appro- 
priate conditions. 



3.4 Generating Asteroid-Sized Masses in Gas 

This section considers the efficient vaporization of significant 
masses of solid matter originally contained in extrasolar as- 
teroids; i.e., processes able to reproduce the rates deduced 
in the previous sections. Such gas production must be orders 
of magnitude more vigorous than sublimation of dust at the 
inner edge of a flat disk (see §3.5), as this cannot produce 
the sufficiently high infall rates dictated by the accretion 
histories represented by Figure [S] The following should be 
considered for heuristic purposes only. 

One possibility is a planetesimal impinging on a pre- 
existing disk fro m the prior asteroid disruptions, as investi- 
gated in detail bv lJural l|2008l ). In that model, if the incoming 
asteroid is less massive than the disk, its entire debris mass 
is reduced to gas via collisions and sputtering, while if it 
is more massive, then the disk mass is vaporized and the 
excess asteroid mass persists in solids. Both cases have the 
potential to produce a large asteroid mass of gaseous debris 
and allow for the possibility that a disk of solids will remain. 
The metal-enriched white dwarfs G166-58 and PG 1225-079 
both appear to have dust rings with enlarged inner holes 



Of tlie expected major elements in a disk of vaporized rock, 
only O I has an ionization potential comparable to H I. 



dFarihi et al]|2010bl . [ioOSbl ). where an impact may have de- 
stroyed the inner region solids but left the outer disk intact. 

A second scenario is a lone tidal disruption event. A 
planetesimal will be shredded on its first approach, but the 
fragments will simply continue along their eccentric orbits 
(i.e., none will impact the star), dispersed but tightly cen- 
tered around the original orbit of the parent body, analo- 
gous to stellar disrup tions at black holes (|Lodato et al.ll200"9l : 
iRosswog et"aLll2009h . The pieces will be sufficiently small as 
to be immune to tidal gravity, but will collide at periastron 
on subsequent passes. Such collisions will occur at high r ela- 
tive velocity and so can generate gas (iTielens et al.|[l993 ) as 
well as dust. Furthermore, solid debris will be prone to sub- 
limation while not yet shadowed by a disk, as optically thin 
material is heated above 1200 K within 1.0 J?© (the approx- 
imate Roche limit) of a typical 15 000 K white dwarf. These 
modes of gas production will eventually abate, as eccentric 
motion damping reduces coUisional velocities and particles 
become shadowed by a growing (flat) disk of material. The 
complete reduction of an asteroid into a disk will likely take 
many orbital periods; e.g., taking the 11.8 yr orbital period 
of Jupiter as a benchmark, one would expect planetesimals 
perturbed from orbits at similar distances from their host 
stars will take hundreds of years to make several passes. 

The timescale for particle collisions in an optically thin 
disk within the Roche limit of a cool white dwarf is typically 
shorter than Poynting-Robertson drag timescale by an or- 
der of magnitude (but with a line ar dependence on optical 
depth and particle size; iFarihi elTal. 2008a ). In this regime, 
planetesimal debris generated by tidal disruption can only 
be reduced by collision or sublimation, at least until the disk 
becomes optically thick, at which point the disk evolution 
will be dominated by Poynting-Robertson drag at the in- 
ner edge, and viscous forces elsewhere (Rafikov. 2011a ). Disk 
observations that reveal 1) silicate emission feat ures in the 
mid-infrared (|jura et al.ll2009l : iReach et al.ll2009l ) and 2) ra- 
dially coincident gas and dust at some polluted white dwarfs 
llFarihi et al. 2012: Molis et al.'2010l: [Brinkworth et"al]|2009l : 
iGansickeet al.. .2008. , .2007 , .2006. ) are consistent with a pic- 
ture whereby the shattered parent body is at least partially 
(if not totally) reduced to micron-sized and smaller parti- 
cles. It is therefore plausible that a significant gaseous mass 
will result from the initial phases of a tidal disruption event. 

3.5 Observational Consequences and Predictions 

This rapid gas accretion scenario has implications for ob- 
servations of circumstellar disks and metal-polluted white 
dwarfs. Perhaps most notably, those stars suspected to be 
accreting gaseous debris alone should exhibit the highest 
rates, and far exceed the Poynting-Robertson drag rates for 
white dwarfs. Ultimately, all material falling onto the stel- 
lar surface will be in gas phase, and circumstellar disks may 
have 1) a cool outer region dominated by solids, 2) a nar- 
row transition region where solids and gas coexist, and 3) 
a hot inner region where all debris is sublimated. Such ar- 
chitecture is inferred fo r the dust disks found and charac- 
terized by Spitzer (e.g ., IFarihi et al1l2009l : I Jura et al.ll2007l : 
Ivon Hippel et aLll2007 , and the limiting factor in their infall 
rates is the (solid) mass per unit time delivered to th e edge of 
the su blimation zone by Poynting-Robertson forces l|Rafikovl 
l2011al ). Once vaporized, the model presented here predicts 
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that gaseous accretion takes place on short timescales, but 
will produce very different infall rates for disks whose total 
mass is dominated by solids, and those dominated by gas. 

Based on the above, those disks detected to have spa- 
tially coincident d ust a n d gas (six are know n at present 
Brinkworth et al.l l201p: iDufour et al.l l2012l: iMelis et al 



20121 : iFarihi et al.ll2012riMelis et al.ll20ld : lBrinkworth et al 



2003 ') are likely to be dominated by solids. Critically, the 



best observations to date indicate that the (calcium) gas 
emission drops off sharply at the inner orbital r adius where 
the particulate disk is model ed to terminate (|Melis et al.l 
l20ld : iBrinkworth et al.1 12009| ). These data are consistent 
with the rapid depletion of a completely gaseous inner re- 
gion, resulting in a relatively low surface density, a corre- 
sponding lack of detected emission, and are compatible with 
the model presented here. 

Notably, observations reveal that white dwarfs with de- 
tectable dust disks have, on average, the highest ongoing ac- 
cretion rates. Figure [4] plots a histogram of metal accretion 
rates for DAZ stars with and without infrared excesses from 
circumstellar disks. While essentially no stars exceed the 
range of infall rates expected for Poynting-Robertson drag, 
the stars without obvious dust are the slowest accretors. 
Many of the DAZ without infrared excess are sufficiently 
warm that they must be accr eting unseen mater ial; either 
gaseous or tenuous dust disks (|Farihi et al.ll2010bh . The pic- 
ture presented here predicts these stars are indeed accret- 
ing from undetected particulate disks, rather than purely 
gaseous matter, but it is also plausib le that a varies substan- 
tially between gas disks (|jurall200^ ). More sensitive infrared 
observations can test this prediction. 

Detailed spectroscopy of DBZ white dwarfs has the po- 
tential to constrain the frequency of events that produce 
significant masses of gaseous debris, or distinct pollution 
events more generally. For a given baseline of exo-chondritic 
or -terrestrial compositions, the heaviest elements such as 
iron and nickel should be under-abundant in DBZ stars if 
the timescale for major asteroid incursions is longer than the 
lifetime of photospheric metals. To date there is a mix of ob- 
servational results: there appear to be two white dwarfs pol- 
luted by intrinsically iro n-poor material (Zuckorm an et akl 
I2OIII : iFarihi et al.ll2011al ) , while most iron-deficient stars are 
consistent with gr avitational settling jKoester et al.l I2OIII : 
IFarihi et al.ll2011bl ). Analysis of a large number of DBZ (and 
very cool DAZ) stars is needed to better understand these 
abundance patterns, as it remains possible that extrasolar, 
solid planetary material can be distinct from that contained 
in any known (or modeled) rocky bodies. But it is plausi- 
ble that certain abundance patterns can only be maintained 
if discrete, substantial pollution events occur on timescales 
shorter than the heavy element lifetimes in typical DBZ 
stars. 

It is noteworthy that the toy model of high and low ac- 
cretion states, which is certainly an approximation of reality, 
underscores the quasi-steady state accretion model recently 
outlined by (|jura fc Xull2012l l. In that model, the heavy ele- 
ment abundances observed in the atmospheres of DBZ-type 
white dwarfs may differ subst antially from both the steady 
state and early phase models ([Koester 2009). Under certain 
conditions, parent body chemical abundances can be uncer- 
tain by a factor of two in mass. The res ults presented her e 
constrain r, but not uj in Equation 10 of ljura fc Xul (l2012h . 



Dust 
Gas? 



5 6 7 B 9 10 11 IS 

log [<dMydt> (g s"')] 

Figure 4. Accretion rates for the 40 DAZ stars in Figure [T] The 
dashed and dotted histograms show stars with and without an 
infrared excess, respectively. Among 19 stars with Tgg > 10 OOOK 
and where ongoing accretion is highly probable, 11 stars accreting 
from dust disks have an average infall rate more than three times 
higher than eight stars that lack infrared excesses. If these latte r 
stars accrete from gas-dominated circumstellar disks l |jurall2008l) , 
the difference in inferred accretion rates would be contrary to the 
predictions made here, which suggest that pure gas disks achieve 
the highest accretion rates, whereas solid-dominate d disks are 
limited by Poynting-Robertson forces l|Rafikovll20"llah . 



3.6 Total Mass Accreted in High State 

An important consequence of the proposed scenario is the 
necessary mass (in gaseous heavy elements) delivered during 
the high state accretion episodes, and specifically the total 
mass accreted within th so that its signature is still visi- 
ble today among the DBZ population. The required mass 
depends on the choice of constants in Equation [S] but for 
r — 1 Myr a mass of Mhth ~ 10^^ g will reproduce the aver- 
age DBZ accretion rates, for Mi equal to the average DAZ 
rate; this solution is represented by solid line in Figure [S] If 
instead the diffusion timescale for the DBZ stars is taken to 
be 10^ yr or smaller (as for GD 362 and Ton 345), then the 
necessary mass decreases to the order 10^^ g. 

If a second planetary body is perturbed to collide with 
a pre-existing disk system, then the gaseous debris mass 
needed to satisfy the above criteria is about that of Vesta, 
conservatively speaking. The largest asteroids in the Solar 
System - Ceres, Vesta, and Pallas - are best described as in- 
tact planetary embryos, and while their sizes and masses are 
exceptional even among the largest objects in the Main Belt 
members, such masses may be more commonplace among 
planetesimals formed at the intermediate mass. A- and F- 
type progenitors of typical, extant white dwarfs. If the entire 
mass is not converted to gas when a second body impacts 
a disk, and remaining material is gradually delivered to the 
star in a particulate disk, then a Vesta analog represents the 
minimum mass necessary to account for the data. 

In the case where singular events are the cause of metal- 
enriched white dwarfs, and substantial gas mass is only pro- 
duced immediately following a tidal disruption, the parent 
bodies must be larger still. Again for the criteria above, one 
needs 10% of a lO^"* g or 1% of a 10^^ g parent body to be 
accreted within a 100 yr period, as gas, with the remaining 
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mass gradually delivered to the star in a particulate disk, 
at the rates observed for the DAZ stars. Here the parent 
bodies would be at least massive as Ceres, and perhaps as 
large as Pluto. Also, the particulate disk lifetimes necessary 
to accrete the remaining mass at more modest rates exceed 
10^ yr. Neither of these predictions is a priori unreasonable, 
and further observations and modeling may help to constrain 
these quantities. 



4 CONCLUSIONS 

The rapid infall of gaseous debris at white dwarfs can ac- 
count for the difference between the currently observed and 
ongoing metal accretion rates at DAZ stars and the time- 
averaged, historical rates inferred for the DBZ stars. The 
data constrain the lifetime of the necessary, high-rate accre- 
tion episodes to be strictly less than 10'* yr, and gas disk 
models with a ~ 0.1 can achieve rates approaching or ex- 
ceeding 10^^ gs~*; sufficient to consume Vesta in 10 yr. The 
several hundred DA white dwarfs surveyed for metals sug- 
gest that the high-rate accretion timescale is likely less than 
10'' yr, while a search of ~ 10 000 DA stars is necessary to 
confidently detect one in a high state that lasts ~ 100 yr. 

Such rapid accretion is not unrealistic and within an 
order of magnitude of the highest accretion rates found for 
cataclysmic variable s in quiescence (|Baskill et al. 2005ll2005l : 
iGansicke et al.iri995l ). Given that ROSAT was capable of de- 
tecting x-ray lum inosities as low as 10^° e rgs"^ at distances 
beyond 100 pc (jPretorius fc Kniggell20l3 ). corresponding to 
accretion rates of 10*^ gs~* for typical 0.6 Mq white dwarfs, 
it is likely that near-future facilities will detect the high- 
rate accreti on episodes predict ed here, if the model is ac- 
curate (see lBear fc Sokerll2012l for alternative predictions). 
eROSITA will be a factor of 100 times more sensitive than 
ROSAT in the 0.5 - 2 keV ban4!|, meaning that accretion 
rates as low as 10**gs~* may be detectable within a few 
hundred pc. Moreover, wide-field and high-cadence surveys 
such as LSST and Pan-Starrs may detect optical transients 
associated with these events. Lastly, GAIA should be capa- 
ble of identifying a single white dwarf accreting in a high 
state, as it should identify 10® white dwarfs and thus yield 
at least 10 000 metal-polluted stars. These data will help to 
constrain th and Mh- 

A readily testable hypothesis is that all relatively warm 
DAZ stars accreting near or below Poynting- Robertson drag 
rates are accretin g from particulate, not gaseous, debris 
disks. iFarihi et all (|2010bl ) has shown that narrow rings of 
dust exist at some polluted white dwarfs, and that 10^^ g 
can easily be missed in photometric surveys for infrared 
excess. Those DAZ stars without obvious infrared excess 
in Spitzer IRAC photometry may reveal dust with higher 
sensitivity observations. High signal-to-noise spectroscopy 
over the 5 — 15 nm range, where both thermal and silicate 
emission are salient for disks observed to date (|Reach et al.l 
I2OO9I : I Jura et al]|2009l ). is perhaps the ideal choice for the 
detection of subtle infrared excesses at white dwarfs. This 
would be somewhat analogous to the identification of dust 
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at HD 69830 (|Lisse et al.l2007l ). where the fractional infrared 
luminosity is relatively low. 

If the accretion of gaseous debris delivers significant 
mass to metal-contaminated white dwarfs, the destroyed and 
devoured parent bodies must be at least as massive as the 
largest Solar System asteroids, and possibly comparable in 
mass to the largest moons. Such objects will be differenti- 
ated, consistent with the few polluted white dwarfs with de- 
tailed abundance measurements for major elements of plan- 
etary solids. The overall results show that quasi-steady state 
accretion may be important, and models that constrain the 
frequency of high accretion states will strengthen the ana- 
lytical connection between the observed stellar abundances 
and the chemistry of planetary debris at DBZ and very cool 
DAZ white dwarfs. 
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